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ABSTRACT 

Using a new grism at the Keck Interferometer, we obtained spectrally dis- 
persed (R ~ 230) interferometric measurements of the Mira star R Vir. These 
data show that the measured radius of the emission varies substantially from 
2.0-2.4 /iiii. Simple models can reproduce these wavelength-dependent variations 
using extended molecular layers, which absorb stellar radiation and re-emit it 
at longer wavelengths. Because we observe spectral regions with and without 
substantial molecular opacity, we determine the stellar photospheric radius, un- 
contaminated by molecular emission. We infer that most of the molecular opacity 
arises at approximately twice the radius of the stellar photosphere. 

Subject headings: instrumentation: interferometers — techniques: interferometric — 
techniques: spectroscopic — stars: AGB and post-AGB — circumstellar matter — 
stars: individual (R Vir) 



1. INTRODUCTION 



Miras are asymptotic giant br anch (AGB) star s that show long-period (hundreds of 



days) photometric variability (e.g.. iFeast et al.lll989l ). These stars typically lose mass at a 
rapid rate (~ 10~ 6 M Q yr~ x ), and this material condenses into heavy elements, molecules, 
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and dust (e.g., IWillsonl 120001 ) . Since AGB stars account for most of the mass loss in the 
Galaxy, their circumstellar environments are the likely factories in which many constituents 
of planets, and life, are produced. 

Separating molecular emission and absorption from stellar photospheric emission is im- 
portant for constraining properties of molecular layers and the underlying stars. The opacity 
and location of gaseous layers constrain the temperature structure in the circumstellar en- 
vironments of Mira stars; this directly affects the dust condensation radius, and in turn the 
mass loss rate. Measurements of stellar radii and temperatures uncontaminated by overlying 
shells of molecular emission or absorption are needed to estimate the location of Mira stars in 
the H-R diagram. These quantities also cons train period -mass-radius relations, with direct 
implications for pulsational properties (e.g., 



Bowen 



1988). 



Historically the circumstellar environments of Mira stars have been studied through high 
resolution infra red spectroscopy, w hich probed layers of CO and H 2 above the stellar pho- 
tospheres (e.g.. iHinkle et al.lll982l ). Interferometric imaging, combined with spectroscopy, 
has more recently allowed spati ally resolved measurements of radio-wavelength co ntinuum 



and spectral line emission (e.g., iReid fc Menten 
dust and molecular emission (e.g., 



Weiner et al 



19971: iDiamond fc Kemballl l2003h . mid-IR 



and molecular emission (e.g., Ilreland et al.ll2004bl ; iThompson et al.ll2002l ). These observa- 



2003) , and optical/near -IR photospheric 



tions showed that molecular emission lies above the stellar photosphere, and that different 
molecules lie at different stellocentric radii. 

Previous near-IR interferometric observations, which probed the hottes t (> 1000 K) cir 



cumstellar emission, used ei ther very low spectral resolution (R < 25; e.g., IThompson et al. 



20021 ; iMennesson et al.ll2002l ) or a few narrow-band fil ters designed to pr obe certain molecular 
lines or nearly line-free regions of the spectrum (e.g.. iPerrin et al.ll2004l ). These observations 
provided evidence for different spatial distributions of stellar continuum and molecular ra- 
diation, but lacked the spectral resolution to clearly separate various components. With 
higher spectral resolution (R > 100), one can distinguish contributions from H2O, CO, and 
the stellar photosphere, and can even discern several of the CO band-heads. Higher spectral 
resolution data therefore allow an estimate of the relative spatial distributions of molecular 
emission and stellar continuum emission in Mira stars. 

Here we present the first observations with a new grism at the Keck Interferometer, 
which enables an order of magnitude higher spectral resolution (R ~ 230 across the K- 
band) than previous studies. We observed the short-period (~ 145 days) Mira star R Vir 
(AKA HP 10994, HR 4808 or IRC 10256). This object varies by ~ 6 magnitudes in th e 



0.6 magnitudes in the K band { Whit cl ock et al 



2000j). 



visible (jKholopov et al. 1985 ) and 
The distance to R Vir is ~ 400-500 pc (lJura &: Kleinmannlll992l ; iKnapp et al.ll2003l ). and 
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the mass loss rate is estimated to be ~ 10 7 M Q yr 1 (I Jural 1 19941 ) 



EXPERIMENTAL SETUP 



We observed R Vir on UT 2006 May 15 using the Keck Interferometer (KI). KI is 
a fringe-tracking lo ng baseline Michelson interferometer combining lig ht from the two 10- 
m Keck apertures (jColavita &: Wizinowichl 120031 ; IColavita et al.l 120031 ) . We measured K- 
band fringes that were spectrally dispersed through a replica grism. The grism has an 
undeviated wavelength of 2.3 /xm in the first order, an apex angle and blaze angle of 6.5°, 
21.36 grooves/mm, and an index of refraction of 1.435 (Infrasil). The KI detector has a pixel 
size of 18.5 /im and a camera with a focal length of ~ 82.5 mm. Thus, the grism provides a 
spectral resolution of R = 230 with 42 10-nm-wide channels across the K-band. 

We measured squared visibilities (V 2 ) for our target and two calibrator stars in each of 
these spectral channels. The calibrator stars are K0-K4 giant stars, with known parallaxes, 
whose K magnitudes are within 0.5 mags of the target. The system visibility (i.e., the 
point source response of the interferometer) was measured with these calibrators, whose 
angular sizes were estimated by fitting blackbodies to literature photometry. The calibrator 
data w ere weighted by th e internal scatter and the temporal and angular proximity to the 



target (IBoden et al 



998h . Source and calibrator data were corrected for detection biases as 



described by lColavital (119991 ) and averaged into 5-s blocks. The calibrated V 2 are the averages 
of 5-s blocks in each integration, with uncertainties given by the quadrature addition of the 
internal scatter and the uncertainty in the system visibility. 

We averaged our data to produce a single measurement of V 2 in each spectral channel. 
Obseravtions of R Vir spanned 0.7 hours, and the averaging has a negligible effect on the 
uv coverage. By comparing our calibrator stars to one another, we see that the channel-to- 
channel variations are small (Figured]). Using HD 107328 (spectral type K0III) to calibrate 
HD 111765 (spectral type K4III), we find a standard deviation of 0.012 in the calibrated V 2 
versus wavelength. Calibrating HD 107328 with respect to HD 111765 we find o = 0.028. 
We therefore adopt a channel-to-channel uncertainty of 2% for our data on R Vir. 

The normalization of V 2 versus wavelength has larger uncertainty than the channel-to- 
channel uncertainties discussed above. Observations of a binary star with a known orbit 
show that the calibrated V 2 have a systematic uncertainty of < 5%. This is a discrepancy in 
the normalization; if we subtract 0.04 from our measurements, the measured and predicted 
V 2 are consistent across all channels at approximately the 1% level. We do not apply this 
offset to our data, but assume that in addition to the 2% channel-to-channel uncertainties 
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described above, the normalization of V is uncertain by ~ 5% 



We used the counts in each channel observed during "foreground integrations" (jColavita 



19991 ) to recover a crude spectrum for our target. We divided the measured flux versus 
wavelength for R Vir by the observed flux from an A1V star (smoothed over the Br7 feature 
at 2.16 /im) and multiplied by a 9000 K blackbody to calibrate the spectral bandpass. While 
our spatially filtered data may contain relatively less flux at shorter wavelengths due to the 
larger effect of atmospheric turbulence, such slopes should be removed through our spectral 
calibration. We estimate that the channel-to-channel uncertainties are ~ 5-10%. However 
there may be larger systematic errors due to different K magnitudes of the target star and 
the A1V spectral calibrator or changing atmospheric or instrumental conditions between 
the source and calibrator observations. We normalized our calibrated photometry and then 
converted to flux units by assuming a K magnitude of 2.0 based on the current visual phase 
and previously studied K-band variability. 



MODELING 



3.1. Star + Shell 



We begin with a simple model consisting of a star surrounded by a single-temperature 
shell whose optical depth varies with wavelength (essentially a Schuster-Schwarzschild model). 
We assume the star emits blackbody emission from a compact continuum-fo rming atmo- 
sphere. We ignore limb darkening, which may affect the stellar radius by ~ 5% (jWelchlll994l ; 



Claret et al. 



19951 ); the effect is comparable to the assumed uncertainties in the data. 



The free parameters of the model are the stellar angular radius and temperature, 9* 
and T*, the shell radius and temperature, # s h e ii and T shell , and the optical depth of the shell 
at each wavelength, t\. Since we must determine the best-fit value of r at each observed 
wavelength, the model has 46 free parameters. We measure V 2 and flux in each spectral 
channel, and we have 84 data points with which to constrain these parameters. 

This model ha s been used previously to model n arrow-band interferometric observations 



of Mira stars (e.g., iPerrin et al.ll2004l ; IWeinerl 120041 ). and can be used to estimate the rela- 



tive spatial distributions of stellar and circumstellar emission. Moreover, comparison of the 
fitted T\ with that ex pected for different molecules (based on laboratory spectroscopy; e.g., 
Rothman et al.l 120051 ) constrains which molecules are present in the circumstellar environ- 
ment. 



We construct the model by considering sight-lines from the center of the star out to 
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the limb of the shell, and summing the contributions from each line-of-sight. The specific 
intensity of the model includes emission from both the star and the shell for sight-lines 

9 < 9*: 



I v (0) = B u (T*)e^/™») + B u (T shell ) [l - e^ 
For 6* < 6 < 9 S h e n, only the shell contributes: 

I u (9) = B u (T shcll ) [i_ e (- 2 Wco S *)] 



(2) 



We compute th e flux and normaliz ed visibility for a series of annuli of infinitesimal angular 
width (see e.g., Eisner et al.l 12004 ) . and sum over all annuli to obtain F u and V 2 for the 
model. We fit our data to the modeled F u and V 2 using a Levenberg-Marquardt non-linear 
least squares fitting algorithm. 



3.2. Star + Molecular Layer 

Instead of considering a shell with arbitrary optical depth t\ (and 42 associated free pa- 
rameters), we can use molecular opacities calculated for stellar atmosphere models (Hauschildt 
private communication). We replace r A with r mo i = N co i x cr mo i(T mol ), where cr mol is the 
molecular cross section (as a function of wavelength) for a single layer in the Hauschildt at- 
mospheric model, and iV co i is the column density of the layer. This column density assumes 
a mix of different molecules (H2, CO, and H2O, among others), and the relative abundances 
are determined through self-consistent modeling of a cool stellar atmosphere. Since the rela- 
tive abundances of H 2 and H 2 and CO may be lower in a shell than in a dense atmosphere, 
our fitted column density for the shell may be over-estimated. The free parameters for this 
model are i?* and T*, the molecular layer radius and temperature, _R mo i and T mo i, and the 
column density of the molecules within the layer, N mo \. 



4. RESULTS AND CAVEATS 

Our modeling indicates that the stellar radius is ~ 1.4 mas (~ 130 R assuming a 
distance of 450 pc) and the stellar temperature is ~ 3700-3800 K (Tabled! Figure [2]). The 
shell of molecular material is cooler and located at a larger radius: -R mo i ~ 2.4 mas (~ 230 
R ) and T mol ~ 1800-1900 K. 

We have assumed in our models that the molecular material around R Vir resides in 
an infinitesimally thin shell. However material is probably distributed over a range of radii, 
and different molecules at different temperatures may contribute to the optical depth. Our 
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simple assumption may be partially responsible for differences between the modeled opacity 
and that expected from molecular line- lists (Figure [3]). Deviations from spherical symmetry, 
not included in our models, may also contribute to differences between models and data. 

The inferred temperature of the molecular shell is substantially higher than the expected 
condensation temperat ure for dust around Mira stars (~ 500 K; e.g.. lRowan- Robinson fc Harris 
19821 ; lHashimotolll994l ) . Dust at larger radii, but within the 50 mas (~ 23 AU) field of view of 
our observations could contribute some over-resolved flux, which would reduce the measured 
V 2 and increase the measured flux. This should be a minor effect given the small amount 
of .fT-band flux from any cool extended dust rel ative to that in the hot, compact molecular 
shell modeled above (see e.g.. Eisner et al.ll2004l ). 



DISCUSSION 



Our analysis suggests that the molecules (predominantly H 2 0) responsible for most of 
the opacity in the circumstellar environment of R Vir form at approximately twice the radius 
of the stellar photosphere. The molecular material cannot be supported by hydrostatic forces 
at this radius. Propagating shocks associated with the pulsation of Mira stars can lead to 
molecular de nsities at these radii orders of ma gnitude larger than predicted by hydrostatic 
models (e.g., iBowenl Il988l ; IWoitke et al.lll999l ). and shocks thus provide a mechanism for 



lifting the molecular layer in R Vir to its observed location. 

The radius and temperature of the molecular shell imply radiative equilibrium with the 
stellar continuum radiation, suggesting that if the shell lies in a post-shock region, it has had 
time to cool. At the molecular densities in such a shell the radiative cooling time is typically 
longer than the time between successive shocks. It is more likely that the molecular material 
in R Vir cooled via expansion; in a rapidly expanding shell associat ed with a hig h mass loss 



rate the cooling time is shorter than the time between shocks (e.g.. lBowenlll988l ). 



Our model fits imply a stellar radius and temperat ure of ~ 130 R^ an d ~ 3800 K. Pre- 



vious broadband measurements of Mira diameters (e.g.. iHaniff et al.lll995l ). which measured 
contributions from stellar photospheres and molecular layers, found large stellar radii, which 
yielded per iod-mass-radius relations that were consistent with overtone mode pulsation (e.g., 
Feastlll996l ). In contrast, the smaller radius determined in this work (and radii determined 
for other Miras from spectrally dispersed or narrow band interferometric observations; e.g., 
Perrin et al.l 120041 ; IWeinerl 12004 ) is consistent with the period-mass-radius relationship ex- 
pected for fundamental mode pulsations, comp atible with other lines of argume nt that Miras 



should pulsate in the fundamental mode (e.g., lBowenlll988l ; IWood et al.lll999l ) 
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A logical next step is to compare spectrally dispersed interferometry data to dynamic, 
self-consistent models of Mira atmospheres that include the e ffects of H9.O opacity. We at- 



temp ted such comparisons but currently-available models (e.g. Jlreland et al.ll2004al ; llreland fe Scholz 



20061 ) have longer periods and are either hotter or have larger photometric variability ampli- 
tudes than our target. These models predict too much flux and not enough H 2 to match 
our data. More appropriate models will become available in the near future (Ireland pri- 
vate communication). Moreover we will obtain additional observations of R Vir, enabling a 
comparison of models and data over a range of pulsation phases 

Data presented herein were obtained at the W.M. Keck Observatory from telescope 
time allocated to NASA through the agency's scientific partnership with Caltech and the 
University of California. The Observatory was made possible by the generous support of 
the W.M. Keck Foundation. We thank the entire KI team for making these observations 
possible, and acknowledge the role of MSC software in our data acquisition and analysis. 
We are indebted to P. Wizinowich, G. Vasisht, and L. Hillenbrand for helping to ensure the 
successful implementation of the grism, and to P. Hauschildt, J. Carr, and M. Ireland for 
providing models and molecular opacities used in our analysis. 
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Table 1. Model Parameter Values 



Model 




0* 


^shell 


#shell 






(K) 


(mas) 


(K) 


(mas) 


(g cm~ 2 ) 


Star+Shell 


3800 ± 440 


1.37 ±0.07 


1800 ± 100 


2.44 ±0.13 




Star+Molecular Layer 


3700 ± 140 


1.40 ±0.07 


1900 ± 100 


2.42 ±0.12 


4.3 ±0.2 



Note. — Uncertainties listed in the table are the quadrature sums of statistical uncer- 
tainties in the fits and assumed 5% systematic errors. 
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Fig. 1. — Calibrated V 2 as a function of wavelength for R Vir and for a K4III reference 
star (note that the reference star is mildly resolved and thus V 2 < 1). The plotted error 
bars correspond to channel-to-channel uncertainties of 2%. We also indicate a 5% systematic 
uncertainty in the overall calibration that affects the normalization (§2D- 
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Fig. 2. — Observed V 2 and fluxes plotted with the predictions of the best-fit models. The 
shaded region shows assumed 10% uncertainties for our flux measurements. Solid and dashed 
lines show the predictions for the star+shell ( §3. 1[) and star+molecular layer ( §3.21) models, 



and contribute substantial opacity from ~ 2.0-2.2 and 
Ro- v ibrational lines of CO contribute opacity between < 



2.3-2.4 /mi (e.g. 



Ludwie 


1971 


"^ 
). 


Rothman et al. 



20051 ). and the four dips in the V 2 and fluxes between ~ 2.3-2 .4 fim occur at the approximate 



wavelengths of the CO overtone bandheads. 
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Fig. 3. — Fitted T\ for the star+shell model (solid line). We also plot r A computed for 
a stellar atmosphere model (Hauschildt private communicati on) in a layer with T = 1800 
K and N m „] = 5 g cm" 2 , and T \ based on line-lists of H 2 (ILudwigl Il97ll ) and CO (from 
HITEMP: [Rothman et aDl2005h with T = 1800 K, A^ H2 o = 5 x 10 20 cm" 2 , and iVco = 10 22 
cm -2 (dashed line). These column densities are not fitted, but rather estimated in an order 
of magnitude sense to produce t\ close to the modeled value. 



